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Abstract 

This article is concerned with an initial boundary value problem for an elliptic- 
parabolic coupled system arising in martensitic phase transition theory of elastically 
deformable solid materials, e.g., steel. This model was proposed in [3], and inves- 
tigated in [3] the existence of weak solutions which are defined in a standard way, 
however the key technique used in [3^ is not applicable to multi-dimensional problem. 
One of the motivations of this study is to solve this multi-dimensional problem, and 
another is to investigate the sharp interface limits. Thus we define weak solutions in 
a way, which is different from [3], by using the notion of viscosity solution. We do 
prove successfully the existence of weak solutions in this sense for one dimensional 
problem, yet the multi-dimensional problem is still open. 

1 Introduction 

In this article we shall investigate an initial-boundary value problem of a new model which 
describes martensitic phase transitions in elastically deformable solid materials, and such 
phase transitions are driven by configurational forces. To formulate this problem, we 
firstly introduce some notations. Let Q be an open bounded domain in with smooth 
boundary dO,. It represents the points of a material body. Define Qt = (0,t) x ^1. We 
use unknown functions: u = u{t, x) is the displacement at time t and position x, T is the 
Cauchy stress tensor, and S is an order parameter which means that if S takes the values 
that are approximately equal to and 1, then the material is in two different phases, say 
7 and 7', respectively. Then the system reads 

-div^r(t,x) = b{t,x), (1.1) 
Tit,x) = D{eiV^u)-eS)it,x), (1.2) 

Stit,x) = -c(-T-e + i>'{S)-i^A,s)\V,S\it,x) (1.3) 
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which must be satisfied in Qt. We prescribe the following Dirichlet boundary and initial 
conditions 



^l[o,t]xan 


= 0, 


(1.4) 


^\[o,t]xdn 


= 0, 


(1.5) 


S\{t=o}xn 


= Sq. 


(1.6) 



In this model, c, u are positive constants, D is the linear, positive definite symmetric 
elasticity tensor. We have chosen the free energy = 7p{e, S, VxS) given by 

V'(e, S, V^S) = ^D{e - eS) ■ (e - eS) + ij{S) + ^| V^^p, (1.7) 

and ^ps is the derivative, with respect to S, of ip. The scalar product of two matrices a, r 
is denoted by cr • r = Yl^j=i ^ij'^ij- There holds the relation 'psi^, S) = —T ■ e + tp'iS). 
e is the strain tensor defined by e = e(Vx.n) = ^ (yxU + ^{^xu)), and the upper-script * 
denotes the transpose of a matrix, e is called the misfit strain. The function ip{S) is 
chosen as a double-well potential for which we assume that 

'ip{S) is smooth and has two minima at S" = and 5 = 1, and 

one maximum at S between and 1, (1-8) 

i)'{S) > 0, if 5 G (0, S) U (1, oo); ^'{S) < 0, if 5 G (S, 1) U (-oo, 0). 

Finally, b = b{t, x) is a given volume force. 

This model was formulated in Q by employing the second law of thermodynamics 
and a formula (see e.g. [H [25} 135) ) of configurational forces. Our model differs from 
the celebrated Allen-Cahn model (which is also called Ginzburg-Landau) by the gradient 
term |Va;5'|. The reason is that in the Allen-Cahn model, the driving force for the motion 
of interfaces is the mean curvature, while the motion of interfaces considered in this paper 
is driven by configurational forces. We mention the key ideas of the derivation. There are 
two main types of phase transition models: sharp interface model and phase field model. 
Our model is derived from a sharp interface model: Assuming that the jump of S, 
across the interface of two phases, becomes smaller and smaller, we see that the equation 
governing the interface approaches to a Hamilton-Jacobi equation St = —c^psl'^xSl which 
is a fully nonlinear equation, thus is difficult to deal with and its solution may develop 
singularities. A usual way for regularizing it is to add an artificial term (for instance, 
lyAxS) as in the theory of conservation laws, but this technique does not work in our 
case. We then think of another type of models, i.e. phase field model, to regularize such 
an equation. To formulate a phase field model, we choose the free energy (|1.7p . and also 
need a suitable fiux which can be chosen in the form 

q = q{ut,T, VxS, St) = T ■ ut + uStV^S. (1.9) 

Then by straightforward computations, we see that if the equations (II. ip - (II. 3p are 
satisfied, then the following Clausius-Duhem inequality is satisfied 

^ V'(e, S, VxS) - dbix q-b-ut<0. (1.10) 
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Hence, we assert the validity of the second law of thermodynamics. For the details of the 
formulation of this model, we refer to the appendix of this paper, or the articles [2l[3ll3]. 

The aim of this article is to propose a suitable concept of weak solutions that works 
for multi-dimensional problem and that makes the investigation of sharp interface limit 
(as V goes to 0) easier, then to prove the existence of such defined weak solutions for 
problem (jl.ip - (jl.6p . There are two most well-known concepts of weak solutions to 
partial differential equations: The first one is the notion of usual weak solutions that 
are defined by employing test functions and the technique of integration by parts, and 
the second one is the conception of viscosity solutions developed by Crandall and Lions 
in 1983, see [15) . etc. In this article, we define a weak solution by combining these two 
notions of weak solutions. To understand why we need two concepts of weak solutions, 
we first investigate the features of this model. Our model consists of a subsystem of linear 
elasticity and a partial differential equation that is degenerate and has strong nonlinearity 
and non-smooth coefficients. The one space dimensional initial-boundary value problem 
for this model has been studied in [3j , in which we define a weak solution in a usual way 
by using a simple technique that makes us possible to rewrite the principle part of the 
equation of the order parameter in a divergence form, i.e. i'Sxx\Sx\ = ^{Sx\Sx\)x- 

However such a technique fails for the corresponding multi-dimensional problem of 
this model, namely vlS.xS\\/ xS\ can't be rewritten in a divergence form. Thus the notion 
of usual weak solutions is not suitable for this problem because we can not reduce the 
order of weak derivatives of the unknown by integration by parts. This is one of the 
difficulties in solving our model. Another one is that the maximum principle, which 
plays a crucial role in the theory of viscosity solutions, is not valid for the whole system 
of equations considered here. So it is not suitable to define weak solutions by using 
the notion of viscosity solutions only. Therefore one of two purposes of this article is 
to propose a suitable notion of weak solutions to this multi-dimensional problem. The 
second purpose is that we shall use our new notion of weak solutions to study, in the 
future, a very interesting problem, i.e. the sharp interface limit of our model. Such a 
problem however may be difficult under the framework of the standard weak solution, 
since the sharp interface problem has a fully nonlinear equation of the order parameter. 

The above consideration leads us to propose a suitable notion of generalized solutions 
to our system by using both notions of weak solutions: we define weak solutions in the 
usual sense for the subsystem of elasticity, and use viscosity solutions to define weak 
solutions to the order parameter equation. Then we construct a sequence of solutions 
to an approximate initial boundary value problem of the system. Applying some com- 
pactness lemma we can show that the limit of the approximate solutions is just weak 
solutions in our sense. Though only the one space dimensional problem is solved up 
to now, we believe this technique works for the multi-dimensional case too. The other 
interesting open problems in this field include: The sharp interface limit of our model, 
and the relationship between weak solutions defined in this article and the ones in 
respectively. 

We are now going to study the definition and existence of weak solutions in a suitable 
sense to problem (jl.ip - (II. 6p in one space dimension, though the definition and some 
a-priori estimates are still valid for multi-dimensional problem. We shall see later on 
that the proof of the existence of weak solutions in this article is significantly simpler 
than that in [3j. 
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Statement of the main result. From now on we assume that all functions only depend 
on the variables xi and t, and, to simplify the notation, denote xi by x. The set = (a, d) 
is a bounded open interval with constants a < d. We write Qt^ := (0,te) x where tg 
is a positive constant, and define 

{v,^p)z = / v{y)ip{y)dy, 
Jz 

for Z = Q oi Z = Qt^. If w is a function defined on Qt^ we denote the mapping 
X I—)- v{t,x) by v{t). If no confusion is possible we sometimes drop the argument t and 
write V = v{t). We still allow that the material points can be displaced in three directions, 
hence u{t,x) S M^, T{t,x) G S"^ and S{t,x) € M, where is the set of 3 x 3 symmetric 
matrices. If we denote the first column of the matrix T{t, x) by T^(t, x) and set 

eK) = ^ (K,0,0) + *K.,0,0)) G 5^ 

then with these definitions the equations (jl.ip - (jl.3p in the case of one space dimension 
can be written in the form 

-Ti = b, (1.11) 
T = Dieiu,)-eS), (1.12) 

5i = c(T-e-4^'{S) + uS^J)\S^\, (1.13) 



which must be satisfied in Qt^. The boundary and initial conditions therefore are 

u{t,x)=0, {t,x) e[0,te] X dn, (1.14) 
S{t, x) = 0, {t, x) € [0, te] X (1.15) 
5(0,x) = So(x), xGJ^. (1.16) 



To define weak solutions to problem (II. lip - ()1.16p , we first introduce some definitions 
on semi-continuous functions. Let / = f{x) be a real function defined in ?7 C M'^ with 
N £ {1, 2, 5}. We denote the so-called upper semi-continuous envelope of / by 

/*(x) : C/^MU{-oo,+cx)} (1.17) 

which is defined by 

f*(x) := limesssup„{/(y) \ y £ U, \x — y\ < r}. (1-18) 

Obviously, f*{x) is upper semi-continuous. And /*(x) := — (— /)*(x) is called lower 
semi-continuous envelope of /. 

We define the Hamiltonian Hx which depends on the unknown T by 

HT{t,x,p,q,r) = c (T{t,x) ■ e- i^'{p) + i^r^ \q\, (1-19) 

where, {t,x) € Qt^, P, Q, r £ M, so (t,x,p,q,r) € M^. It is easy to show that if T is a 
continuous function in {t, x) and tjj' is continuous in S, then we have that Ht is continuous 
in {t, x,p, q, r), thus 

{HTy{t,x,p,q,r) = {HT)*{t,x,p,q,r) = HT{t,x,p,q,r). (1.20) 
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We now can introduce the notion of weak solutions for our problem. In what follows 
we shall assume that p is a real number such that 

p>l. (1.21) 



Definition 1.1 A function {u,T,S) which satisfies that 

{u,T,S) G L~(0,te;^o(^^)) X L'=^{0,t,;L\n)) x L°"(QtJ, (1.22) 
is called a weak solution to system \1.11\) - il.l6\) if 

I) for almost every t E [0,te], equations (jl.lip . (jl.l2p and the boundary condition 
(I1.14P are satisfied weakly. 

II) S is a viscosity solution to equation \1.13\l . if S satisfies both i) and ii) below: 

i) S is a sub-viscosity solution to equation \1.13\) . i.e. for any function (p{t,x) 
in C^'^(Qte), if S* — (p attains its local maximum at {T,y), then 

MT,y) < {HT)4'r,y,S*{T,y),(j)^{T,y),(j),^^{T,y)), (1.23) 

and S*{0,x) < So{x); 

ii) S is a super-viscosity solution to Eq. f fi.3|) . i.e. for any function 4>{t,x) in 
C^'^(Ote); if — (j) attains its local minimum at {T,y), then 

Mr,y) > {HTnT,y,S,{T,y),Mr,y)AxxiT,y)), (1.24) 

and S^:{0,x) > So{x). 

Now we are able to state our main result as follows. 

Theorem 1.2 Suppose that b, bt G C{[0,te]; L^{il)) for any given positive constant te, 
and that Sq € Hq{^). Furthermore, we assume that the function satisfies the assump- 
tion lll.8\} . 

Then there exists a weak solution (n, T, S) to problem - il.l6\) in the sense of 

Definition \l.l\ and in addition to ^.22\) . we have that the solution satisfies 



Our notion of generalized solutions is a combination of the concept of usual weak 
solutions and the notion of viscosity solutions. This idea comes partly from some discus- 
sions with Prof. Alber and partly from the paper by Giga, Goto and Ishii |24j which is 
concerned with the global existence of weak solutions, however without uniqueness, to 
the system consisting of a semi- linear diffusion equation in two disjoint open sub-domains 
denoted by ^±{t) of one simply connected domain J7 (The complement of union of these 
two parts is so-called the interface r(t)), and a nonlinear interface equation. The system 
is composed of the interface equation 

V = W{v)-CK, on r(t) (1.25) 
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and the diffusion equations 

vt = i^Au + g±{v), for x E ^^±(t), t > 0. 

Here, V = V{t,x) is the speed of T(t) at a; € r(t) in the normal direction of n from 
to ^}-{t). K is the mean curvature of T{t) at x G r(t), v is the density. And 
c, are positive constants, W,g± are given bounded continuous functions over R. Note 
that in the work |24) , the driving force for the motion of an interface is due to the mean 
curvature (see formula (I1.25P ). while the motion of an interface considered in this article 
is driven by configurational forces and the motion is governed by V^[«S'] = cn ■ [E]n (the 
sharp interface case), where E is the Eshelby tensor, an energy- momentum tensor, see 
[H pp. 753-767]. 

We recall the literature related to our results. There have been many papers on the 
theory of viscosity solutions since the notion of viscosity solution was proposed in 1983 
by Crandall and Lions [_15j. This notion is applicable to the scalar equations or the 
weakly coupled systems, for which the maximum principle holds. Hence, the comparison 
theorem is valid, this plays an important role in the proof of uniqueness of viscosity 
solution. For an overview of the theory, we refer for instance to Capuzzo Dolcetta and 
Lions [12J, Ishii and Lions ^29j, Crandall, Ishii and Lions ^4j, Jensen[31j, Crandall and 
Lions ^6j, Ishii [27j, Souganidis [H] for the scalar equation case, and to Engler and 
Lenhart |18] . Ishii and Koike [28], etc. for the system case, and the references are cited 
therein. For the background of our model and mathematical results related this article, 
we refer the reader to work by Alber and/or Zhu [21 O HI [5l [H [71 |l2], Kawashima and 
Zhu [32], Ou and Zhu [37]. 

The main difficulties and our strategies in the proof of Theorem 11.21 are as follows: 
Firstly, the definition of weak solutions is a new problem since our system comprises of 
a linear elliptic system of u and a nonlinear equation of S which can not be rewritten 
in the divergence form. Secondly, the equation for the order parameter is degenerate 
and its coefficients is not smooth. To overcome these difficulties, we make a suitable 
smooth approximation of the non-smooth term which leads the equation of the order 
parameter to a uniformly parabolic equation with smooth coefficients. We employ the 
energy estimates to discuss the limits of approximate solutions. 

The remaining of this article is organized as follows. In Section 2 we state an ap- 
proximate initial boundary value problem, and apply the existence theorem in the book 
by Ladyzenskaya et al. |33] to prove existence of classical solution to this approximate 
problem. Then we derive in Section 3 the uniform a priori estimates which are indepen- 
dent of a small parameter k for the approximate solutions. Then we apply the a priori 
estimates, a lemma of the Aubin-Lions type and a theorem on the stability of viscosity 
solutions to discuss the limits and prove the existence of weak solutions in the sense of 
Definition 1.1. Finally Section 4 we present briefly in the appendix the derivation of our 
model. 

2 Existence of solutions to the modified problem 

In this section, we are going to study an approximate initial-boundary value problem and 
show that it has a classical solution for any fixed positive constant k. Since we shall let 
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K go to zero, we may assume, without loss of generality, that 

< K < 1. 

Let X ^ Cq°(M?, [0, oo)) be a function satisfying xi^i x)dtdx = 1. We set 

, . 1 ( t X 
Xnit^x) := —X ( -, - 

and for h e L°°(Qt^,M) we define 



{XK*b){t,x)= XK{t - s,x - y)b{s,y)dsdy. (2.1) 
Jo 

We smooth the term 15^1 as follows 



\S.U = V\WT^, (2.2) 
and choose a sequence such that 

SS^C^m, \\SS-So\\HHn)^0 (2.3) 

— >■ since Cf^{Q) is dense in Hq{Q). 
Then the smoothed initial boundary value problem of (II. lip - (11.16P turns out to be 

-ri = x.*b, (2.4) 

T = D{e{u,)-eS), (2.5) 

St = cu\Sx\kSxx + c(t ■ e - ip'{S)) {\Sj:\f, - k). (2.6) 



and the boundary and initial conditions become 

u\[o,t,]xdn = 0, (2.7) 
S\[o,te]xdn = 0, (2.8) 
•S'l{o}xQ = ^0- (2.9) 

By the choice of Sq, we see that the compatibility condition |on = is met. 



Remark 2.1. There are some other ways, which are different from (|2.2p . to smooth 
the function \p\. We need only to require that the smoothed equation (12. 6p for the order 
parameter meets the assumptions of the maximum principle. 

To prove the existence of classical solution to the approximate problem (12. 4p - (12. 9j) . 
we employ the Leray-Schauder fixed-point theorem (see, e.g. [33]) to the following prob- 
lem 

-ri = Xx.*b, (2.10) 
T = D{e{ux)-XeS), (2.11) 

St = cu\SxUSxx + c(T-e-i;'{S)^{\Sx\n-K). (2.12) 
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Here A € [0, 1]. The boundary and initial conditions are 



u\[o,te]xdn = 0, (2.13) 
5|[o,te]xac = 0, (2.14) 
•S'l{o}xn = ^Sq. (2-15) 

Define for any S € S := C^+t'2+"(QtJ (here < a < 1) a mapping Pa : [0, 1] x i3 ^ 
B; S 1-^ S where 5 is a solution obtained by the following procedure: 

i) For any fixed S, it is easy to find a unique solution (u, T) which depends on S, to 
the following boundary value problem for almost every given t 

-Tl = Xx.*b, (2.16) 
T = D(e{u,)-XeS^ , (2.17) 
u\dn = 0. (2.18) 

ii) Then inserting this T into equation (I2.12|) we can obtain a unique classical solution 
S to problem (l2T2]l . and (l2T5]l . 

With the help of some a priori estimates, we then obtain 

Theorem 2.1 Suppose that all the assumptions in Theorem \1.2\ are met, and the com- 
patibility conditions Sq = Sq^x = 5'o,xx = at x = a,d are satisfied. 

Then for any fixed K > 0, there exists a unique classical solution {u,T,S) G C'^'^{Qt^)x 
C^'i(QtJ X C2+"'i+"/2(QtJ to problem ^KM " WM) which satisfies 

St. G L\QtJ. (2.19) 

Remark 2.2. The compatibility conditions in Theorem 2.1 are different from usual ones 
and they are derived as follows: From the system and initial data, there must hold 

r(0,x)U=,,,- Z)e(n.(0,x))|,^,,, = 0, 

(2.20) 

v\Sq^x\kSq^xx + T{Q,x) ■ £{\Sq^x\k - K')\x=a,d = ^■ 

Note that the values ofux{0,x) at boundary can be arbitrary, so isT{0,x). Thus from the 
definition of the function | • 1^ we see that the second term of ^2. 20\) is satisfied provided 
that Sq^x = Sq^xx = at x = a,d. 

The proof of Theorem 12.11 is equivalent to that P\ has a fixed point for A = 1. It is 
clear that if A = 0, Pa maps any S € B to S = 0, i.e. PqS = 0. Thus by the Leray- 
Schauder fixed-point theorem, we see it remains to derive the estimates that are stated 
in Lemma 2.2 - Lemma 2.5 and Lemma 2.7, from which we conclude the compactness of 
Pa. To obtain those a priori estimates, we assume that there exists a classical solution 
{u,T,S) G C^^HQtJ X C^'HQtJ X C2+°'i+"/2(QtJ to problem (l2T0]l - (f2T5]l such that 
Stx € L^iQtJ. 

Firstly, writing T in terms of S we reduce system (I2.10p - (12.12P into a scalar equation 
with a non-local term, and apply the maximum principle this equation to obtain 
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Lemma 2.2 There holds for te > 

In this lemma and in the fohows context, we denote by C a constant which is in- 
dependent of K, but may depend on v, while a constant C may depend on both k and 

Proof. To make use of the maximum principle, we solve {u, T) in terms of S from the 
first two equations (j2.10p - ()2.1ip . provided that S is given. Then the whole system can 
be reduced into a single equation, but with a nonlocal term. We need some notations as 
used in |4j. Let be the subspace of all matrices A ^ with Aij = for i,j = 2, 3. The 
orthogonal space to <S^ is denoted by . It consists of ^ G <S^ satisfying An = An = 
for all i = 1,2, 3. Note that £{ux) € <S^. Let P be the canonical projection of 5^ into 5^. 
Since D : 5^ ^ 5^ is a positive definite linear mapping, {a, r) = Da ■ r defines a scalar 
product on 5^. The projection of onto 5^, which is orthogonal with respect to this 
scalar product is denoted by Q. These definitions imply that 



Define further that 



we then obtain 



kerQ = D-'^S^ = D-\crP. 
e* = Qe, u* = (eJi,2e*i,2£*J, 



uit,x) = Xu*(l S{t,y)dy - j-^ I S{t,y)dy ] + w{t,x), (2.22) 



d 



De* 

T{t,x) = D{e* -e)XS{t,x) - / XS(t,y)dy + (j(t,x), (2.23) 

d-a Ja 

where the function {w{t,-),a{t,-)) (here t is temporarily regarded as a parameter) is the 
unique solution of the following boundary value problem 

-o-i^(x) = b{x) in n, (2.24) 
a{x) = De{w^{x)) in n, (2.25) 
w{a) = f{a), w{d) = f{d) (2.26) 

and b = Xb{t), / = 0. Note that u* € M^, e* E depend only on the misfit strain e. 
Inserting the formula of u, T into equation (j2.12p yields that system (|2.10p - (j2.12p is 
reduced into a single equation for S with a nonlocal term. Invoking the definition of 
we see that the assumptions required by the maximum principle are satisfied. Thus we 
can apply the maximum principle to this single equation and the proof of this lemma is 
complete. 

Next we can derive the following estimates for the derivatives of S. 
Lemma 2.3 There holds for any t € [0, te] that 

\\SAt)f + f [ \SxUSa:.fdxdT < C, (2.27) 
Jo Jn 

j^[i\SxUS,,\)l + \St\l) dxdr < C. (2.28) 



'0 Jn 

2 



Here and hereafter, we denote the L -norm over O by 
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Proof. By definition we have the property \p\f^ > k, from which we obtain 

0< \pU-k< \pU + k<2\pU. 

Using estimate (j2.2ip and formula (|2.23p . recalling the assumptions on 6, one concludes 
that 

\\T\\l^(q,^) < C. (2.29) 
Note that Stx G L'^iQt^), for any fixed k, implies that 

2^11*^2:11^ = {Sx,Sxt)- 

Multiplying (I2.12|) by —Sxx and integrating the resulting equation with respect to x, 
using integration by parts, and invoking the estimates (j2.2ip and (|2.29p we get 

^■^WSxW"^ + ci^ {\Sx\kSxx,Sxx) = c(^(t ■ e - ip'{S)^ {\Sx\k -!<■),- Sx 

< C'^|5'x|I, I "Sa:; 1 1 1 iSxx I ^ ) (2.30) 

where we used the notation {f,g) = f{x)g{x)dx. Applying the Cauchy-Schwarz in- 
equality, we infer from ()2.30p that 

Id - i i 

tj;itII'S'x|| + cu {\Sx\kSxx,Sxx) < C\\\Sx\K\\L^{n)\\\Sx\ K Sxx 1 1 

< C{\\Sx\\'^ +l)\\\Sx\ISxx\\ (2.31) 

By the Young inequality and the property that < jp| + k for any k > 0, we derive 
from ([OT]) that 

2^il'5'x|| ~^ CU {\Sx\kSxxt Sxx) ^ i I 1 5*3; 1 K 5'za:: 1 1 + C'i^( || 5*2; || + 1) 



< ^ 

- 2 



/ \SxU\Sxx\^dx + C4Sxf + C. (2.32) 



Thus we arrive at 



j^WSxf + CV j^\SxUSxx?dx < C^Sxf + C. (2.33) 

Using the Gronwall inequality to (j2.33p one can easily obtain (|2.27p . 

By the interpolation technique and (|2.27|) . we have that for some 2>p>l,q = | 
and 1 + ^ = 1 that 

t f 

V 



/o 




Jn 



< 



< 



{\Sx\K\Sxx\f dxdr 
{\Sx\.)^ ({\Sx\S^ \Sxxf) dxdT 

{ f I {\Sx\S'^ dxdrV ( f I {\Sx\S'^ \SxxrdxdT\ 
\Jo Jn J \Jo Jn / 

[ [ ilSxU^dxdr) ' ( [ [ ISxUSxxl^dxdrY . (2.34) 
Jo Jn J \Jo Jn / 



' t 
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Invoking the property that \p\^ < \p\ + k and inequahty (|2.27|) yield that for < 2, 



i.e. p < i, the right hand side of (j2.34p is bounded. 



Making use of (|2.34|) (with p = I) and equation (|2.12p we have for any test function 



— ^111 I K I 



i'ISxIkSxx + (T • e - ip'{S)){\S,j, 



+c 



\Sx\\ + i)|Iv'|Il4(q^j 



— ^ ( II I'S'xLiS'a:; 



'LT!{Qt, 



(2.35) 



< C\\ip\\L4(^Q^^), 

where we apphed the Holder and Young inequalities. Thus we arrive at \\St 
thus prove (|2.28p . And the proof of this lemma is complete. 

For the solution to the elliptic part of the system, i.e. ()2.16p - ()2.18p . we have 
Lemma 2.4 There hold for almost every t G [0, tg] that 

Mt)\\w^,p^n) + \\Tit)\\LP(n) < C, (2.36) 



uit)\\H^n) + \\TmHHn)<C. 



(2.37) 



Proof. Using the estimate (I2.2ip . we get S{t) G L^(il) for almost every t G [0,te] since 
the domain is bounded. Recalling estimate ()2.27p . we obtain easily (|2.36p - ()2.37p . 
by the regularity theory of elliptic systems (or just using the formula (I2.23j) since our 
problem is one dimensional). This completes the proof of the lemma. 

Now we differentiate (j2.10p once formally with respect to t and recall the assumption 
on bt, then use again the theory of the elliptic system to get 



Lemma 2.5 There hold for almost every t G [0, t^] that 



'-t\\LP{n) 



<C(l + \\St\ 



mn)), ll^tlli§(Q,j 



T eC{[0,te];C'^in)), and ||r| 



< c, 

< c. 



(2.38) 
(2.39) 



To prove the above lemma, we shall make use of the following lemma which is of 
Aubin-Lions type, see, for instance, Lions ^34J, and for the case r = 1, see Simon [39] , 
Roubicek 



Lemma 2.6 Let Bq, B, Bi be Banach spaces which satisfy that Bq, Bi are reflexive 
and that 

BqCCB C Bi. 



Here, by CC we denote the compact imbedding. Define 

df 



W 



f \feL°-{0,t,;Bo 



dt 



eL'{0,t,;Bi: 



with te being a given positive number and 1 < r < oo. 
Then the embedding of W in C{[0,te\; B) is compact. 
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Proof of Lemma 2. 5. We need only to prove ()2.39p . From (12.380 and Lemma 12.41 we 
assert that 

T e L^{0,te;H\n)), and Tt e lI {0,te; lI {Q)). (2.40) 
Thus we can choose 

Bo = H\n), i3 = C"(0), Bi = Ll{n), r = ^, 

which meet the requirements of Lemma 12.61 and a € (0, Whence (I2.39P holds. And 
the proof of Lemma 12.51 is complete. 

Furthermore, for any fixed k, we have 
Lemma 2.7 There hold for any t G [0,te] that 



\\St{t)f + 




Jn 



\Sx\n + '^)\Sxt\'^dxdT < C, 

\\s,At)\\ < c. 



(2.41) 
(2.42) 



Proof. We prove firstly ()2.4ip . To this end, we differentiate equation ()2.12p formally with 
respect to t, then multiply the resulting equation by St and integrate it with respect to 
X to get 



't\?-cu {\Sx\.Sxx)tStdx-c ({T-e-ij'{S)){\SxU-K)] 5* dx = 0.(2.43) 



ld_ 

2dt' 

It is easy to see that 



I I K SxX ) I 



\i\K.di ) Sxtdx 



\i\Kdi 



xt 



\Sx\K\Sxt\ dx. 



(2.44) 
(2.45) 



Thus ()2.43p becomes 



+ CV \SxU\Sxt?dx - c ({T-e-^'{S)){\SxU-n)] Stdx = 0. (2.46) 



1 d 
2(ft 

We now handle the last term on the left-hand side of (j2.46p as 
T-e-i^'{S)^ {\SxU-'i))^Stdx 
Tfe-4'"{S)St {\Sx\^ + l)\St\dx + C [ T-e-i,'{S) \{\SxU)' SxtSAdx 




< c 



< cmw + \\st\\) (II \SxUL^in) + mstw + c\\st\\ wsxtW 

< C( (1 + II ISxUL^^n)) \\Stf + (II \Sx\.\\L^(n) + 1)11^*11 \\St\\) + ^\\Sxtf. (2.47) 



Here we used the estimates ||(T, 5)||loo(q^j < C, \\Tt\\ < C(l + \\St\\), |(|?/U)'| < C, and 
the Holder, Young inequalities. By the Sobolev imbedding theorem, we have 



II \SxitML^in) < Ci\\Sx{t)\\Lo.in) + 1) < Ci\\Sx{t)\\H^n) + !)• 



(2.48) 
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Applying the estimate p.38p . which is vahd for p = 2, combining ()2.47p and (I2.46|) . we 
arrive at 

j^\\St(.t)f <C{1 + + C + l) . (2.49) 

We shah make use of the Gronwall inequahty in the fohowing form 

y'{t)<A{t)y{t) + B(t) imphes y(t) < y(0)e^' ^^"^^ + / B{s)e^s Mr)dT (2.50) 

Jo 



where y, A, B are functions satisfying that y{t) > 0, A{t)^ B{t) are integrable over [0, tg]- 
Defining 

L\ iio/iMi2 , no II \ ^ /'i , II ii2 



y{t) = \\St{t)\\\ A(t)=C(l+||5,|lHi(n)), S(t) = C 1 + 115,11^1 



from the estimate ||S'a;||/^i(Qj^-) < C which is a consequence of (|2.27p and the fact \p\^ > k, 
it fohows that the above defined A{t), B{t) are integrable over [0, tj- Thus we can apply 
(|230]1 to ([TiQ]) and obtain 

\\St{t)f<C, 

whence 

\\St{t)f + f I \S,\^\S^t?dxdT < C. (2.51) 
Jo Jn 

From which we obtain easily (j2.4ip since k at this moment is a given number. Therefore 
we can use the equation to get easily (I2.42p . Thus the proof of this lemma is complete. 

Remark 2.3. Since we use the imbedding (j2.48p which is valid only in one dimensional 
case, thus this lemma is only true for this one dimensional problem. 

Remark 2.4. To derive (j2.4ip rigorously, we employ the technique of finite difference as, 
e.g., in J77| ^. We assume that there exists a unique classical solution {u,T,S) to problem 
(|2T0]1 - (|2T5]1 such that 

{u,T,S) E C''\QtJ X C''\Qt^) X C2+"'i+-/2(Q,J, S^t G i'(Qte). 

Define Sh{t,x) = {S{t + h, x) — S{t,x))/h for any h>0. Then from (I2.12p we obtain 



ru f rSx(t+h,x) \ ^ (t+h,x) 

Sm = -T-i \^Ud^) +r{T-e-i,'{S)){\S,\^-K) , (2.52) 

^ \Js4t,x) I ^ h (t,x) 



for any [t € [0, t^^ — 5], where 5 is a fixed number such that 5 > h. Here and hereafter, we 
use the notations = f{t + h,x) — f{t,x) and f\(^t,x) = f{t}^)- Multiplying ()2.52p 

by Sh and integrating the resulting equation with respect to t, x over Qt yield 

fte—S r 1 f-Sx(t+h,x) 

\\Sh{t)f+ci^ / / - / iCl^d^Shxdxdt 



n ^ JSxit,x) 



|2 



\\Sh{0)r + c / / {Th-e-r{S*)Sh){\SxU-'^)\it+h,x)Shdxdt 
Jo Jn 
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Here S* is a number between S{t + h,x) and S{t,x). Note that the second term on the 
left hand side of (|2.53p is equal to 



cv 



I / \Sx{Vix)\^\Shx?dxdt>C\\ShAh(Qt^_iy 



where rj £ [t,t + h] and we used \p\i^ > k. By definition, we have 

\Sxit + h,x) + Sx{t,x) 



\Sx{t + h,x)\^ + \Sx{t,x)\ 



< 1. 



So the second and third terms on the right hand side of ()2.53p are of lower orders and 
can be estimated in a similar way to (j2.47p . We thus arrive at 

Further, we write 

l-te-S 







te-S 



\Sx{r],x)\^\Shx\'^dxdt 

((|5,(r/, x)U - \Sx{t, x)U) + \Sx{t, x)U) \Shx?dxdt. (2.54) 



Invoking the Holder continuity of Sx, applying the Fatou lemma for any fixed 5, we take 
the limit as h ^ Then letting 5 — )• 0, we justify (j2.4ip and omit details. 

We now turn back to prove Theorem 12.11 
Proof of Theorem \2.1\ To complete the proof of the global existence of classical solution, 
we need to prove that ||5a;|Lc«/2,Q((5, ) < C. To this end we make use of the estimates 
listed in Lemmas 12.21 - I2.5l and Lemma 2.7. To prove this, we invoke the following lemma 
see, e.g. [33] 



Lemma 2.8 Let f{t,x) be a function on Qt^ such that 

i) f is uniformly (with respect to x) Holder continuous in t, with exponent 

< a < 1, that is \ f{t,x) - f{s,x)\ < C|t-s|", 

and 

a) fx is uniformly (with respect to t) Holder continuous in x, with exponent 
< /3 < 1, that is \ fx{t,x) - fx{t,y)\ < C'\y - x\P . 

Then fx is uniformly Holder continuous in t with exponent 7 = a(3/{l + (3), such that 

\fx{t,x)-fx{s,x)\ <c"\t-s\'', yxen,o<s<t<te. 

where C" is a constant which may depend on C,C' and a,f3. 

By applying this lemma we assert that there exists a constant < a < 1 such that 
||5a;||cc«/2,c, < C. By the a priori estimate of the Schauder type for parabolic equations, 
we thus obtain that 

\\S\\(jl + a/2,2 + a{Qt^) < C, 

which ensures us to apply the Leray-Schauder fixed point theorem, and the proof of global 
existence of classical solution is complete. Using the technique of difference quotient with 
respect to t to this classical solution see e.g. [33] we can prove ()2.19p . And the proof of 
Theorem 12.11 is thus complete. 
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3 Existence of weak solutions 



3.1 Uniform a priori estimates 

This subsection is devoted to derivation of some uniform a priori estimates, which are 
independent of n € (0, 1], for the approximate solution to ()2.10p - (j2.15p . However these 
estimates may depend on z/, this thus makes it difficult to discuss the sharp interface 
limit J/ — >■ 0. To investigate such a sharp interface limit, we need new techniques. 

We now denote the approximate solution by {u'^, T'^, S'^). Therefore we collect a priori 
estimates, which have been established in Section 2 and are independent of k. 

Lemma 3.1 There hold for any t G [0,te] that 

\\SS{t)f + f [ \S^U\SS.\'dxdT < C, (3.1) 
Jo Jn 

l^({\S^US^JP + \Snl)dxdT < C, (3.2) 

|2 

Remark 3.1. From (|2.3ip we see that the constant C depends on v. 



'0 Jn 

ft 



S^Wmmdr < C. (3.3) 







3.2 Limits 

With the help of Lemma 2.6, applying the uniform a priori estimates established in 
Subsection 3.1, we shall investigate in this section the limits, as k — > 0, of the approximate 
solutions and complete the proof of Theorem 11.21 

Proof of Theorem \1.2l Firstly we apply again Lemma 12.61 to show that the sequence of 
the approximate solution S'^ has a subsequence which converges strongly. To this end, 
we choose 

and 

1 4 
< a < -, Pi = 3- 

It is easy to see that such defined Bq, Bi are reflexive. Therefore we apply Lemma 12.61 
and conclude that the sequence {S'^}^ is a compact in C{[0,te];C°'{Q)). Thus we can 
select a subsequence of it, and denote it by {5"^"}„, such that, as n ^ 00, 

Kn 0, 

and 

5"" ^5, inC([0,te];C°(n)), (3.4) 

from which we obtain that 

115"" - 5'||c([o,t^]xn) 0- (3-5) 
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On the other hand, by Lemma 2.5, we assert that there exists a subsequence of T"^ such 
that 

^ T in C7'^(gtJ, (3.6) 
from this, (j3.5p and ()2.22p . we obtain consequently that 

{u^^,T^^) (u,T), uniformly in C'^COtJ, (3.7) 

as n ^ oo. 

In what follows, we are going to prove that the limit function (n,T, S) is just a weak 
solution to problem (jl.lip - (jl.l6p in the sense of Definition 11.11 It is not difficult to 
show that (ll.lip and (I1.12P are satisfied by the linearity of those two equations and the 
uniform convergence of u'^, T'^. The remaining part of this section is to prove that (jl.l3p 
is satisfied. 

We shall make use of the theorem on the stability of viscosity solutions, see e.g. [23j. 

Theorem 3.2 (Stability of viscosity solutions) Assume that Fn converges to F lo- 
cally uniformly (as n — )■ oo) in the domain of definition of F. Assume that Vn is a 
viscosity solution to 

{vn)t + Fn{t,x,Vnyvn,V'^Vn) < {resp. > 0) in Qt^, 

and that w„ converges to v locally uniformly in Qt^ as n ^ oo. Then v is a viscosity 
solution to 

vt + F{t, X, V, Vu, V'^v) < [resp. > 0) in Qt^ . 
To apply this theorem, we define Vn = S'^" and 
Fnit,x,p,q,r) = HT'^n{t,x,p,q,r) = c (T''"{t,x) -e-Tp'ip) + i^r^ {\q\^„ - k„). 

Invoking (j3.7p and (j3.5p we conclude that 

i) S'^" converges to S locally uniformly in any compact subset in Qt^. 

ii) HT^-n (t, x,p, q, r) converges to Ht locally uniformly in any compact subset in (0, te) x 
X M X M X M. 

iii) Since S""" is a classical solution to equation (12.12P when T'^" is regarded temporarily 
fixed, S'^" is also a viscosity solution to ()2.12p . 

Therefore, we can apply Theorem 4.1 and conclude that the limit is a viscosity solu- 
tion to St = Hxit, X, Sx, Sxx)- Hence, recalling the properties S*{t, x) = ^^(t, x) = S{t, x) 
and {HT)*{t,x,p,q,r) = {HT)*{t,x,p,q,r) = HT{t,x,p,q,r), we assert that {u,T,S) is 
a weak solution to problem (jl.lip - ()1.16p in the sense of Definition 11.11 and the proof 
of Theorem 11.21 is thus complete. 
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4 Appendix 



Since our model is quite new, we briefly sketch, for the sake of readers' convenience, the 
physical background and the derivation of the diffusive interface model (jl.ip - (II. 6p from 
a sharp interface model. We also refer the reader to [31 HI [5]. Our model differs from the 
Allen-Cahn model by a gradient term. The main reason is: In the Allen-Cahn model, 
the driving force for the motion of interface is the mean curvature, while the motion of 
interface considered in this article is driven by configurational forces, see e.g. [25| [35]. 

Material phases are characterized by the structure of the crystal lattice, in which the 
atoms are arranged. An interface between different material phases moves if the crystal 
lattice in front of the interface is transformed from one structure to the other. Often phase 
transformations are triggered by diffusion processes. A well-known model for diffusion 
dominated transformations is the Allen-Cahn equation when the order parameter is not 
conserved (or the Cahn-Hilliard equation if the order parameter is conserved) . We derive 
our model (II. ip - ()1.6p from a sharp interface model for diffusionless transformations, also 
called martensitic transformations, see e.g. [26', p. 162]. This sharp interface model is 
an initial-boundary value problem for the unknown functions u, T and for the unknown 
interface T{t) (ZVt between two material phases, which is a free boundary. It consists of 
(jl.ip - (II. 2p and the interface conditions 

V{t,x)[S]{t,x) = c(^-{T){t,x)■e[S]it,x) + [^l,{S)]{t,x)^ , (4.1) 

[u]{t,x) = 0, [T]{t,x)n{t,x) = 0, (4.2) 

which must hold for x G r(i), and of a Dirichlet boundary condition for u the initial 
condition (jl.6p . We use the notation [/] = f+ — f~ and (/) = |(/+ + /-), where /+, /_ 
are the limit values of the function / on both sides of T{t). Moreover, V(t,x) € 
denotes the normal speed of the interface T(t), which is measured as positive in the 
direction for which [S]{t,x) is positive. Here c is a positive constant. Equation (14. ip . a 
constitutive equation, determines the normal speed V of the phase interface as a function 
of the term — (T) • e[S] + [ip{S)]. Some computations show that this term is equal to the 
expression n-[E]n with the Eshelby tensor E (an energy- momentum tensor, see |19. p753- 
p767]) and the normal vector n to T{t) (cf. 0]) and thus is a configurational force. We 
assume that V depends linearly on the configurational force, which is the most simple 
constitutive assumption. Thus, in this model the evolution of the phase interface is 
driven by the configurational force along the interface, an assumption appropriate for 
martensitic transformations. 

Though configurational forces were introduced in the first half of the last century, 
it was clearly stated for the first time in [1] that (jl.ip . ()1.2p . ()4.ip . ()4.2p form a closed 
initial-boundary value problem. Applications of this model can be found, for example, 
in im [36| HO] , where equilibrium configurations for materials with phase transitions are 
determined, and in [30] , where the evolution of phase interfaces in ferroelectric materials 
is modeled. In a sense, this free initial-boundary value problem from solid mechanics is 
comparable to the Stefan problem in fluid mechanics. 

The initial-boundary value problem (jl.ip - ()1.6p can be considered to be a regulariza- 
tion of this sharp interface model, which could be used to prove existence of solutions of 
the sharp interface model, and it can also be considered to be a diffusive interface model 
for martensitic phase transitions, which is useful by itself and avoids some disadvantages 
of the model with sharp interfaces. We are interested in both aspects. 



17 



The derivation of (II. ip - ()1.6p given in [21 0] uses a rigorous method. To make the 
model plausible, we derive the model here in a different, short, but formal way. To this 
end we replace the phase interface T(t), across which the order parameter jumps from to 
1, by finitely many interfaces parallel to the original interface, and consider a new order 
parameter, again denoted by S, with small jumps across these interfaces, such that the 
sum of the jumps is equal to 1. We assume that the new order parameter satisfies (|4.ip 
and (14. 2p along all interfaces. If we increase the number of interfaces and decrease the 
jump height, the new order parameter will converge to a continuous or even differentiable 
order parameter, for which the normal speed of the level manifolds is equal to the limit 
of the normal speed of the interfaces. For this limit speed we obtain from (14. ip 

V{t,x) = clim{{T)-e + i;'{S*)) = c{-T-e + i^'{S)) = cMe{V^u),S). (4.3) 

The limit order parameter thus satisfies the Hamilton-Jacobi transport equation 

St = -ciPs{e{S/^u),S)\V,S\, (4.4) 

since the level manifolds of solutions of equation ()4.4p have this normal speed. The 
idea suggests itself to approximate the solution of the sharp interface model by smooth 
solutions {u,T,S) of the system (jl.ip . (II. 3p . (14. 4p . Yet, examples in one space dimension 
show that in general the function S in such a smooth solution develops a jump after finite 
time. The reason for this is that the function ip' appearing in ips is not monotone, since 
V' is a double well potential. After S has developed a jump, ()4.4p can no longer be used to 
govern the evolution of S. To avoid this problem and to force solutions to stay smooth, 
(j4.4p has been replaced by (jl.Sp . which contains the regularizing term i^l'S/ xS\AxS with 
the small positive parameter z^. This yields the model p.ip ~ (II. 6p . 

The choice of this special regularizing term follows from the second law of thermody- 
namics, which every model must satisfy. This law requires that there exist a free energy 
Ip and a fiux q such that -^ip + div^^ q < b ■ ut holds; cf. [8j. If we choose a free-energy 
and a flux as (jl.7p and ()1.9p . it follows by a short computation for solutions {u,T,S) of 
(fTTI) . (fOl) that 

d 

—Ip - divx {Tut + i^StV^S) - b-ut = {ipsie, S) - i/A^S')S't. 
ot 

Inserting (jl.6p into this equation shows that the right-hand side is non-positive, whence 
the second law is fulfilled. However this would not be true by using, as in the theory of 
conservation laws, the standard regularization (i.e. adding an artificial viscosity term) 
St = -cM^'^xu), S)\V^S\ + lyA^S of 
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